Faithful Effective- One- Body waveforms of 
small-mass-ratio coalescing black-hole binaries 
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We address the problem of constructing high-accuracy, faithful analytic waveforms describing the 
gravitational wave signal emitted by inspiralling and coalescing binary black holes. We work within 
the Effective-One-Body (EOB) framework and propose a methodology for improving the current 
(waveform) implementations of this framework based on understanding, element by element, the 
physics behind each feature of the waveform, and on systematically comparing various EOB-based 
waveforms with "exact" waveforms obtained by numerical relativity approaches. The present paper 
focuses on small-mass-ratio non-spinning binary systems, which can be conveniently studied by 
Regge-Wheeler-Zerilli-type methods. Our results include: (i) a resummed, 3PN-accurate description 
of the inspiral waveform, (ii) a better description of radiation reaction during the plunge, (iii) a 
refined analytic expression for the plunge waveform, (iv) an improved treatment of the matching 
between the plunge and ring-down waveforms. This improved implementation of the EOB approach 
allows us to construct complete analytic waveforms which exhibit a remarkable agreement with the 
"exact" ones in modulus, frequency and phase. In particular, the analytic and numerical waveforms 
stay in phase, during the whole process, within ±1.1% of a cycle. We expect that the extension of 
our methodology to the comparable-mass case will be able to generate comparably accurate analytic 
waveforms of direct use for the ground-based network of interferometric detectors of gravitational 
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I. INTRODUCTION 



Coalescing black hole binaries are among the 
most promising gravitational wave (GW) sources for 
the currently operating ground-based detectors like 
GEO/LIGO/ VIRGO. The most useful part of the wave- 
form for detection comes from the most relativistic part 
of the dynamics: the last few cycles of the adiabatic inspi- 
ral, the plunge and the merger. Since LIGO is currently 
taking data at the expected sensitivity, it becomes ur- 
gent to have accurate template waveforms for detection. 
In the most general, spinning case, these waveforms are 
complicated functions of the initial masses mi and mi 
and the spins Si and S2 of the two constituent black 
holes. Due to the multi-dimensionality of the parame- 
ter space, it seems impossible for state-of-the-art numer- 
ical simulations to densely sample this parameter space. 
This motivates a renewed effort for developing analytical 
methods towards computing the huge banks of accurate 
template waveforms needed to densely cover the param- 
eter space. 

As far as we know, the first estimate of the complete 
waveform (covering inspiral, plunge, merger and ring- 
down) of a coalescing, non-spinning black hole binary 
was made in 2000 on the basis of a new analytical 
approach to the general relativistic two-body dynamics, 
the Effective-One-Body (EOB) approach (2j. At the time, 
no reliable numerical simulations of coalescing black hole 
binaries were yet available. This EOB-based estimate 
used a 2.5 post-Newtonian (PN) accurate description 



of the dynamics (the only one available at the time 1 ) 
and a quadrupole-type "restricted waveform" approxi- 
mation for the waveform down to a "matching radius" 
''match where the two-black-hole system was replaced by 
a unique, ringing Kerr black hole. It was also argued that 
the matching radius r matc h could be taken as the EOB- 
deformed "light ring", tlr, which was noticed to be very 
close to the point where the orbital frequency reached a 
maximum. As a first approximation, a complete wave- 
form was computed by matching, at r matc h — ?"lr, the 
inspiral-plus-plunge waveform to a ring-down waveform 
made of the least damped quasi-normal mode (QNM) of 
a Kerr black hole of mass and angular momentum corre- 
sponding to those of the relative dynamics at r matc h- In 
other words, Ref. [l| was making the assumption that the 
merger phase of the two black holes was very brief, and 
did not correspond to any especially strong GW emission 
feature, but would simply smoothly connect the plunge 
behavior to the ringdown one. 

In 2001, a combination of (rather short) full numeri- 
cal simulations with a "close-limit" approximation Q led 
to the first numerical estimate of the waveform emitted 
by the plunge from the ISCO, followed by merger and 
ringdown [ 5j, [6( . Very recent breakthroughs in numeri- 
cal relativity simulations have finally made accessible a 



1 This EOB waveform estimate was later updated by taking into 
account the more accurate 3.5 PN dynamics, as well as by con- 



sidering spinning black holes Q|. 
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much more precise knowledge of the waveforms generated 
by the merger of two black holes of comparable masses 
mi and 7712, possibly with spin. These results have been 
obtained by different groups with different approaches, 
and exhibit convincing internal con verg ence and a nice 
mutual consistency [J 1, ©, EE ED, H El El El El- 
The availability of reliable numerical results makes it ur- 
gent to compare analytical and numerical approaches for 
(i) understanding in depth the physics that is involved in 
the process and thus (ii) completing the currently avail- 
able analytical knowledge of the problem with the "miss- 
ing" non-perturbative physics provided by numerical rel- 
ativity simulations. Several different ways of combining 
the knowledge acquired through analytical and numer- 
ical methods have been recently explored, such as the 
use of "hybrid" signals made by stitching together an in- 
spiral PN-type waveform to a NR merger one [13, Ell , or 
the s tudy of whether some simple multiparameter (BCV- 
type [19() analytical waveform can have sufficiently large 
overlaps with hybrid signals. Our aim here is different. 
We wish, ultimately, to construct high-accuracy purely 
analytical template waveforms, depending only on phys- 
ical parameters, and covering the full process: inspiral, 
plunge, merger and ringdown; by high-accuracy we mean 
that their phasing (and possibly their amplitude too) 
would be very close to the true phasing, uniformly, to 
a few percent. Let us recall that it is convenient (fol- 
lowing the terminology of Ref . [20( ) to distinguish effec- 
tual and faithful templates: effectual templates are de- 
fined by the property of having sufficiently large (say 
> 96.5%) overlaps with the expected "real" signals, af- 
ter maximization over all (kinematical, dynamical, and 
possibly fudge) parameters. Effectual templates maybe 
useful tools for detecting real signals, but they might 
lead to large biases in the measurement of the physi- 
cal parameters of the system. By contrast, faithful tem- 
plates are defined by the property of having not only 
large overlaps, but also of being so "close" to the real sig- 
nals that the overlap between the template and the real 
signal is maximized for values of the parameters which 
are very close to the real ones ("small biases"). Recent 
"first order" comparisons between full numerical relativ- 
ity waveforms and EOB-based waveforms have already 
established that EOB waveforms can be effectual (over- 
laps > 96.5%) over the full mass range of interest for 
GEO/LIGO/ VIRGO, say a total mass M = m Y + m 2 
between 10 and 120 solar masses [13, [H|. Our main 
aim will be to prove that EOB-based waveforms can also 
be faithful over the full mass range of interest. Our 
methodology for achieving this aim is to improve the 
current waveform implementations [3, H, [I?], HH of the 
EOB approach 0, H, [H, [H H3] by understanding, el- 
ement by element, the physics behind each feature of 
the waveform, and by comparing EOB-based waveforms 
with "exact" waveforms obtained by numerical relativity 
approaches. In the present paper, we shall tap the in- 
formation contained in small-mass-ratio waveforms. [A 

preliminary version of our results was reported in [25j.] 



Indeed, the limiting case mim^ <C (mi+1112) 2 constitutes 
a "clean laboratory" where many subtle physical issues 
can be studied in detail (and by means of rather light nu- 
merical tools), without the complications entailed by the 
current three-dimensional numerical simulations (resid- 
ual eccentricity, numerical noise,...). We rely on the re- 
cently developed numerical waveform computations a la 
Regge- Wheeler- Zerilli in Ref. [26j. Note that we consider 
non-spinning inspiralling systems circularized by radia- 
tion reaction. 

The various improvements in the EOB waveform con- 
struction that we shall address here concern: (i) an im- 
proved analytical expression for the ((£,m) = (2,2)) 
waveform which includes a resummation of the tail ef- 
fects, and a 3PN-accurate "non-linear" amplitude cor- 
rection, (ii) the inclusion of non-quasi-circular correc- 
tions to the waveform, (iii) a comparative study of the 
two radiation reaction expressions during the plunge that 
have been proposed ([l| versus [24|), (iv) the inclusion of 
non-quasi-circular corrections to radiation reaction, and 
(v) an improved treatment of the matching between the 
plunge and ring-down waveforms (which takes into ac- 
count a new understanding of the importance of the num- 
ber of QNMs, the sign of their frequencies, and the length 
of the interval on which the matching is done) . We shall 
show that the resulting improved implementation of the 
EOB approach yields very faithful waveforms whose am- 
plitude and phase agree remarkably well with the "ex- 
act" ones: in particular, the EOB phasing differs from 
the "exact" one by less than ±1.1% of a cycle over the 
whole process. 

This paper is organized as follows. In Sec. |TT] we re- 
call the EOB construction for the relative dynamics of 
the binary system and specify it to the small mass ratio 
limit. In Sec. [ill] we discuss the computation of the GW 
emission and we focus on the comparison between differ- 
ent expressions of radiation reaction. Sec. IIVI introduces 
our "best" EOB-type waveform for the transition from 
inspiral to plunge and shows that it is an excellent ap- 
proximation to the "actual" waveform. Some conclusions 
are presented in Sec. [V] 



II. EOB RELATIVE DYNAMICS 

In this section we recall the EOB non-perturbative 
construction of the two-body dynamics including a ra- 
diation reaction force and then specify it to the small 
mass limit. The EOB approach to the general rela- 
tivistic two-body dynamics is a non-perturbatively re- 
summed analytic technique which has been developed in 
Refs. Q,i,[l|22|,[2l,[I3. This technique uses, as basic in- 
put, the results of PN theory, such as: (i) PN-expanded 
equations of motion for two point-like bodies, (ii) PN- 
expanded radiative multipole moments, and (iii) PN- 
expanded energy and angular momentum fluxes at infin- 
ity. For the moment, the most accurate such results are 
the 3PN conservative dynamics [13, HI, and the 3.5PN 
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energy flux [29j, |30|, |31| . Then the EOB approach "pack- 
ages" this PN-expanded information in special resummed 
forms which extend the validity of the PN results beyond 
the expected wcak-field-slow-velocity regime into (part 
of) the strong-field-fast-motion regime. At the practical 
level, the result of this "packaging" is that the compli- 
cated PN-expanded relative dynamics, (in the center of 
mass frame) of a binary system of masses toi and mi 
is mapped into the simpler (essentially geodesic) dynam- 
ics of a particle of mass /j = m\mij{m,\ + m 2 ) moving 
in some effective background geometry (in Schwarzschild 
gauge) 

ds 2 = -A(r)dt 2 + B(r)dr 2 + r 2 (d9 2 + sin 2 0dip 2 ) . (1) 

Here and below we work with dimensionless reduced 
variables r = R/M and t = T/M, with M — m\ -t- 
mi\ (r,9,ip) are polar coordinates in the effective prob- 
lem that describe the relative motion. The dynamics 
of the effective particle is described by a Hamiltonian 
H~eob(M, /i) and a radiation reaction force J-"eob(-M, t 1 )- 
In the general comparable-mass case -Heob has the struc- 



ture H E ob(M,h) = M\ll + 2v{H v -l) where v = 

n/M = m\mij{m\ + m 2 ) 2 is the symmetric mass ratio. 
In the test mass limit that we are considering, v « 1, 
we can expand -ffEOB m powers of v. After subtract- 
ing inessential constants we get a Hamiltonian per unit 
(/x) mass H = lim^o^EOB — const. )//i = \im u ^o H u . 
It is convenient to replace the Schwarzschild radial co- 
ordinate R by the Regge- Wheeler tortoise coordinate 
i?„ = R + 2M\og[R/{2M) - 1], and, correspondingly, 
the radial momentum Pr by the conjugate momentum 
Pn t of R* . Then the (specific) Hamiltonian reads [26[ 



H= x lA(r)[l + ^)+pl 



(2) 



Here we have introduced dimensionless variables r = 
R/M, r* = R*/M, p rt = Pr Jh and Pip = P f /(/j,M), In 
the v — > limit, we have that A(r) = i3(r) _1 = 1 - 2/r. 
We have used the tortoise canonical pair (r* , p r , ) instead 
of the "traditional" one (r,p r ) in the Hamiltonian for 
two related reasons. On the one hand, p r ^ has a fi- 
nite limit when r tends to the zero of A(r) (r = 2, the 
event horizon), while p r diverges there 2 . On the other 
hand, since one of our goal is to compute the gravita- 
tional wave signal from the relative dynamics, it is very 
useful to have the (radial-delta-function) source moving 
quasi-uniformly towards — oo on the doubly-infinite r* 
axis, rather than exponentially slowing down and getting 



stuck at r = 2M (this allows us to cut-off more smoothly 
the radial motion after some time). For details of our 
numerical implementation see [26[. 

Hamilton's canonical equations for (r, r* , p rr ,p v ) in the 
equatorial plane (9 — n/2) yield 
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Note that the quantity Q, is dimensionless and represents 
the orbital frequency in units of 1/M. In these equa- 
tions the extra term (of order 0(u)) represents the 
non conservative part of the dynamics, namely the radi- 
ation reaction force. During the quasi- circular inspiral, a 
rather accurate expression for J- v is the following Pade- 
resummed form 



(8) 



which is expressed in terms of the PN ordering parameter 
Vq = il 1 / 3 . In this expression, the function /dis denotes 
the "factored flux function" of [2(| (scaled to the New- 
tonian (quadrupole) flux), taken here in the v — ► limit. 
Ref. [l| (in the comparable mass case) assumed that the 
analytical continuation of the expression (jHJ) might still 
be a sufficiently accurate description of radiation reaction 
effects during the plunge. On the other hand, the authors 
of Ref. [24J pointed out that the expression (J8J) assumed 
the continued validity of the usual Kepler law f2 2 r 3 ~ 1 
during the plunge. [This is why we label the expression 
([5]) with a superscript K, for Kepler.] They, however, 
emphasized that the Kepler combination K = fl 2 r 3 sig- 
nificantly deviates from one after the crossing of the Last 
Stable (circular) Orbit (LSO), to become of order of 0.5 
at the light ring. Ref. [24[ went on to argue for a differ- 
ent expression for the radiation reaction, say T v (without 
any superscript) that does not assume Kepler's law. In 
the v — > limit this new expression reads 



f 

T = f 

J~ in — 

A* 



5 1 - y/2n v 



(9) 



2 This divergence leads, in particular, to numerical problems in the 
source of the Regge- Wheeler-Zerilli equation: some "exact" ana- 
lytical cancellations that are explicit when expressing the source 
in terms of p r * (which stays finite), become implicit when us- 
ing p r and numerically tend to yield a problematically divergent 



where the ordering PN parameter is given by the az- 
imuthal velocity v v — rQ,. Note that the essential dif- 
ference between the two possible expressions for the ra- 
diation reaction is that cx f2 7 / 3 , while T v cx Q 5 r 4 . 
For simplicity, in both possible expressions for the radi- 
ation reaction our current estimate of /dis is obtained 
by Pade approximating the 2.5 PN results, a procedure 
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FIG. 1: Relative dynamics for v = ji/M — 0.01 and initial 
separation r — 7. The plot shows the (nearly indistinguish- 
able) trajectories with (black line) and with T v (red line). 
The dashed circle indicates the LSO (which is crossed at the 
dynamical time t/2 ~ 240). The light ring (at r = 3), near 
which the plunge waveform will be matched to a ring-down 
one, is indicated (dotted circle) . It is crossed at the dynamical 
time t/2 ~ 300 J2(J. 



III. REGGE-WHEELER-ZERILLI WAVEFORMS 

Following what we did in [261 ] . the computation, a la 
Regge- Wheeler- Zerilli 3 of the gravitational waveform in 
the v — > limit needs three separate steps: (i) to initialize 
the system dU)-®, (ii) to integrate it in time and (hi) to 
solve the Regge- Wheeler and Zerilli-Moncrief equations 
as an initial value problem, with source terms driven by 
the particle dynamics. In practice, once the dynamics is 
computed, one needs to solve (for each multipole (£, m) 
of even (e) or odd (o) type) the following couple of de- 
coupled partial differential equations 

f)2\V,( e /°) ~ 2 ^(e/o) tr(e/o)^(e/o) _ q(e/o) ^ 

with source terms Sj^J^ linked to the dynamics of the bi- 
nary. A thorough discussion of the analytical and numer- 
ical details of our approach has been given in Ref. [26| . 
so that it is not necessary to repeat it here. We only 
recall that the initial condition for the relative dynam- 
ics is given as in (26| . notably by specifying a non-zero 
initial value for p r> ( "post-adiabatic" approximation). In 

addition, the gauge-invariant functions are propor- 

tional to the actual GW polarization amplitudes, and so 
their knowledge encodes the complete information about 
energy, momentum and angular momentum losses. 



that has already been shown to be quite robust versus 
"exact" numerical results in 20]. 

Let us finally comment on the numerical values of the 
parameter v that we shall consider. First, we should 
consider values such that the ^—dependent corrections 
in the two-body effective Hamiltonian H Vl i.e. the 
i'— dependent contributions in the EOB metric functions 
such as A(r) = 1 — 2u + 2vu z + ■ ■ ■ (where u = 1/r), are 
indeed numerically negligible. However, we do not wish 
to consider too small values of v. Indeed, the present 
study of the inspiral, plunge and merger of a mass \i 
into a Schwarzschild black hole of mass M is intended as 
a toy model of the coalescence of two comparable-mass 
black holes. We therefore wish to minimize the quanti- 
tative difference between these two types of inspiralling 
motions. In particular, it has been shown in Ref. [l| that 
the number of orbits during the comparable-mass plunge 
was roughly ~ (4^) -1 / 5 , i.e. a number of order one when 
mi = rri2 (i.e. v = 1/4). We should not therefore con- 
sider values of v so small as to make (4i/) -1 / 5 much larger 
than one. As a compromise we shall always consider here 
the value v = 0.01, which gives only ~ (4z/) -1 / 5 ~ 1.9 
orbits during the plunge (as will indeed appear in Fig. [1] 
below) . 



A. Angular momentum loss 

Let us consider the relative dynamics of a particle of 
mass ix — 0.01M plunging from r = 7. The relative or- 
bit (for both choices of radiation reaction) is shown in 
Fig. m The dashed circle indicates the LSO, which is 
crossed at t/2 ~ 240, while the dotted circle indicates 
the light ring (LR), which is crossed at t/2 ~ 300 [26| . 
The first point we want to discuss is the comparison (and 
contrast) between the two expressions of the angular mo- 
mentum loss mentioned above, the "usual" Keplerian one 
of Eq. (J5J and T v of Eq. © . These two expressions 
are equivalent during the inspiral due to the validity of 
Kepler's law f2 2 r 3 = 1 along circular orbits, but they 
start strongly deviating from each other after the cross- 
ing of the LSO. We wish to know how they compare to 
some sort of "exact" angular momentum flux during the 
plunge, and whether one of them is to be preferred. 

Let us now be precise in which respect we can compute 
an "exact" result for the angular momentum loss during 
the plunge. We recall that via Eq. (17) of Ref. [26[ one 
can compute the instantaneous angular momentum flux 
at infinity dJ/dt (as a sum of multipoles) from ^llf 1 - 
We have done this for both T/^ and T v including all 



3 i.e., by computing the metric perturbations of a Schwarzschild 
spacetime by means of a multipolar decomposition. See for ex- 
ample [32ll and references therein. 
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the multipoles of the radiation up to £ — 4. The solid 
lines in Fig. [5] refer to this computation. There, dJ/dt is 
extracted at r obs = 500 and is shown versus observer's 
retarded time u = M(t° hs — r° hs ) + Au, where one shifts 
the retarded time u by a constant so as to best fit (dur- 
ing the early inspiral) the dynamical time t which enters 
as the argument of the local angular momentum loss of 
the source, i.e. the radiation reaction T^(t) or T^it). 
We use Au = 0.4M. The labelling "Flux with cx ft 7 / 3 
Rad Reac" indicates the outcome of a dynamics with 
(black line), while "Flux with cx fl 5 r 4 Rad reac" 

refers to a dynamics using T v (red line). Fig. [2] shows 
that the differences in the fluxes at infinity are really 
tiny. In that sense, we know with good accuracy (assum- 
ing balance between losses at infinity and local losses) 
what should be the "exact" (local) angular momentum 
loss of the system. By contrast, the two dashed lines in 
Fig. [3 represent the two analytically-assumed local an- 
gular momentum losses; i.e., the two radiation reaction 
terms T^(t) or P v (t). The good news is that they agree 
quite well with the "exact" result during not only the in- 
spiral, but also a fair fraction of the plunge (which starts 
at u/2M ~ 240). However, the bad news is that they 
start differing significantly from the exact result (as well 
as from each other) during the later part of the plunge 
(say, after u/2M ~ 290, remembering that the matching 
to ring-down will be done around u/2M ~ 300, which 
corresponds to light-ring crossing). The "usual" Kepler- 
based radiation reaction T~ is too large (as anticipated 
in Ref. [24| which pointed out that K — ft 2 r 3 sharply 
decreases during the plunge, if we notice that the ratio 
fiy jTtp ~ K~ 4 / 3 ), while the hopefully better radiation 

reaction T v turns out to be somewhat too small towards 
the end of the plunge. On the other hand, the near coin- 
cidence of the two corresponding fluxes at infinity shows 
that the "analytical uncertainty" in the radiation reac- 
tion does not matter much for computing the waveform. 
This confirms, and makes more quantitative, the finding 
of [H, [33[ that, after the crossing of the LSO, the relative 
motion is essentially quasi-geodesic; i.e., relatively little 
influenced by the radiation-reaction force. 

The physics behind the differences between the "ana- 
lytical" and "exact" angular momentum fluxes (as well 
as between the two analytical ones) lies in the presence 
of non- quasi- circular (NQC) effects during the plunge 
phase. During the adiabatic inspiral one takes advantage 
of the quasi-circularity (QC) of the motion to operate 
two separate simplifications: (a) one neglects the terms 
proportional to dr/dt , and (b) one neglects the terms 
proportional to d 2 r/dt 2 , which, by using the equations 
of motion, give (modulo the square of dr/dt) terms pro- 
portional either to v 2 — l/r ( "virial combination" ) , or to 
f2 2 r 3 — 1 ("Kepler combination" minus 1). During the 
plunge, both simplifications become less and less valid 
so that one should correct any QC-derived expression by 
terms proportional to powers of the two basic dimension- 
less NQC ratios (dr / dt) / (fir) and (l/r - v 2 )/v 2 . When 



Flux with Q " Rad Reac 

- Flux with r 4 Rad Reac 

- - Rad Reac « £1 7 ' 3 

- - Rad Reac « £1 5 r" 

Rad Reac <* r 4 with non-QC correction 




50 100 150 200 250 300 350 

U/(2M) 

FIG. 2: Comparison between various angular momentum 
losses: GW fluxes at infinity (solid lines; computed a la Regge- 
Wheeler-Zerilli including up to £ — 4 radiation multipoles, 
using two types of radiation reaction in the driving dynam- 
ics: 'K' (black) or 'not K' (red)); or mechanical angular mo- 
mentum losses: (dJ/dt) K — —P^ ( hyphenated black lines), 
dJ/dt — —Tip (hyphenated red lines). The dotted line refers 
to the mechanical losses with non-quasi-circular corrections. 
See text for discussion. 



considering a time-even NQC correcting factor, it should 
have the general form 

where we used the same "Kepler combination" K = ft 2 r 3 
as before. In our present case, we physically expect the 
presence of such a time-even NQC 4 factor in the radia- 
tion reaction. As we expect (from the usual quadrupole 
formula, as well as from the well-known fact that or- 
bital eccentricity tends to increase the flux of angular 
momentum) that the NQC factor in the radiation reac- 
tion should be larger than one, we conclude that the "im- 
proved" radiation reaction expression T v proposed in [24j 
is a physically more appealing starting point for applying 
such a NQC correcting factor. In addition, as, pragmati- 
cally speaking, the two types of terms parametrized by a 
and a' in Eq. above behave in a roughly similar way 
during the plunge (they are both positive and they both 
increase), we can try, in first approximation, to use only 
one of them, say the first one parametrized by a. 

In Fig. [5] we show (as a black dotted line labelled "Rad 
Reac oc f2 5 r 4 with non-QC corrections") the evolution of 



4 Note that we can henceforth view NQC as meaning either non- 
quasi- circular or next-to-quasi-circular. 
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the NQC-corrected radiation reaction 



jrNQC 



T 9 \ 1 



RR Pr, 



(12) 



with the numerical value a RR — +2.9 (fixed by eye). We 
see how such a simple (and physically reasonable, because 
larger than one) NQC factor suffices to have a good agree- 
ment between a "NQC-corrected analytical radiation re- 
action" and the "exact" angular momentum flux dJ/dt 
up until the light-ring, that is crossed at u/(2M) w 300. 
This gives our first instance of the aim of the present 
paper: to better understand (by comparing to "exact" 
numerical results) the "missing physics" in some first- 
order analytical expression, and then to use this improved 
understanding to further improve the starting analytical 
expression. Note that, for any starting QC expression X 
that one wishes to correct, it is in principle possible to 
compute in General Relativity the "real" NQC correct- 
ing factor fx^ C , of the form (fTTj) with some analytically 
computable coefficients ax, a' x , etc. Note, however, that 
one expects these coefficients to become (when consider- 
ing PN corrections) functions of 1/r and other dynamical 
variables. As in several other known cases, one expects 
that these PN corrections might significantly affect the 
"effective" numerical value of the coefficients ax, a' x , 
etc needed for a good representation of the considered 
quantity during the plunge, ft would be interesting to 
have analytic (at least lPN-accurate) results on the co- 
efficients a RR , a' RR to compare them with our presently 
"experimental" estimate a R g^ ctivc = +2.9. 



B. Exact waveforms 

Let us now consider the waveforms obtained by nu- 
merically integrating the Regge- Wheeler- Zerilli equation 
(TTU1) . with a source corresponding to the plunging dy- 
namics discussed above, with the radiation reaction 
(we have checked that the NQC-improved radiation reac- 
tion (jlip introduces only very minor changes in the wave- 
form). For the most relevant multipolarity [I = m = 2, 
with even parity), the resulting "exact" waveforms (de- 



noted ^ 



(exact) 



or sometimes ^22) are displayed in Fig. [3] 
The left panel of the figure exhibits (black lines) the 

real (top) and imaginary (bottom) part of '3> 2 ' 2 xact '' (nor- 
malized as in [32j), The waveforms are extracted at 
fobs = 500 and are shown versus the (shifted) observer 
retarded time. They exhibit the usual "chirp-like" struc- 
ture corresponding to the inspiral, and end up with the 
black hole QNM ringing. 

The most useful information about the dynamics of 
the merger hidden in the (exact or analytical) waveform 
can be extracted by inspection of the corresponding in- 
stantaneous gravitational wave frequency. Since ^22 is a 



complex number, this quantity is computed as 



M<jJ 2 2 = W 2 2 = -9 



*9 



(13) 



where 3 indicates the imaginary part and the overdot 
stands for d/dt. The exact u> 22 is shown as a black solid 
line in the right panel of Fig. [3] In the same plot we 
also include, as a blue line, twice the orbital frequency O 
(which is the zeroth-order analytical prediction coming 
from any "standard quadrupole approximation"). This 
plot shows that ufa***^ ~ 2fl during the inspiral and 
most of the plunge (the LSO is crossed at u/ (2M) ~ 240). 
It is only when the (relative) radius gets quite near the 
light-ring value r = 3 (at u/(2M) 300, indicated 
by a green dashed vertical line) that the two frequen- 
cies start deviating from each other. While the dou- 
bled orbital frequency 2Q reaches a maximum equal to 
2f2max = 0.2703 very near the light ring 5 and then expo- 
nentially decreases, the exact GW frequency i2> 2 2 keeps 
growing after the light ring until it saturates and os- 
cillates around the fundamental QNM frequency of the 
black hole, Mu 220 = 0.37 3 67 1 5 [Mill. The oscillation 
in the GW frequency around the fundamental QNM fre- 
quency is (essentially) due to the superposition of both 
positive- 6 and negative-frequency excitations in the ring- 
down waveform. 7 Let us mention that this kind of be- 
havior of MjfJ^ (odd and even-parity GW frequencies) 
is absolutely general and the amplitude of the oscillation 
depends on m [36J . 

The fact that the amplitude of the frequency oscilla- 
tion in the right panel of Fig. [3] is quite small then means 
that the positive-frequency mode is excited with a much 
larger coefficient than the negative-frequency one. In 
view of our general aim of understanding all the sepa- 
rate physical elements crucially entering the final wave- 
form, it is interesting to explain what is the basic under- 
lying physics behind this large dissymmetry in the exci- 
tation coefficients of the positive- frequency and negative- 
frequency modes (we have checked that this dissymmetry 
extends to the first few overtones above the fundamen- 
tal QNM mode). Though this basic physics can be de- 



5 It is easy to see from the Hamiltonian equation J5} why this is so: 
neglecting the relatively slow variations of the angular momen- 
tum and the energy, Q vary proportionally to the ratio A(r)/r 2 
which is precisely the effective radial potential for photon (or 
ultra-relativistic) circular orbits. 

6 As in Quantum Field Theory we call positive frequency a mode 
oc e -1 "' with ui > 0. We also use for the QNMs (see below) 
the notation oc e~ CTt . In this notation, the positive frequency 
fundamental mode has Oq~ = c?o + i^O 1 where ojq indicates the 
frequency and 00 the damping time of the fundamental mode, 
while the negative frequency one has = ag — iujq. 

7 Evidently, had we initialized the dynamics with the opposite sign 
of p v , i.e., the particle inspiralling clockwise and not counter- 
clockwise, we would have obtained a mirror-reflected result, with 
the negative frequency modes being more excited than the posi- 
tive frequency ones. 
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FIG. 3: Left panel: Quadrupole waveforms (t = m = 2, even parity): exact Zerilli-type (black line) and "matched" improved 
EOB-type (red line). Right panel: exact and matched-EOB modulus and instantaneous gravitational wave frequency. The blue 
line displays twice the orbital frequency Q. The two vertical dash-dot lines mark the ends (at 297.8 and 301.4) of our matching 
interval (which is centered, at 299.6 on the maximum of the modulus, and which includes the light ring, at 300.9). 



scribed in terms of equations, it is best understood by 
means of some diagrams. One should have in mind that 
solving the Regge-Wheeler-Zerilli equation (fTOj) means 
essentially considering a point-source (moving along the 
radial axis) whose "strength" oscillates proportionally to 
exp[— imtp(t)] (with m = +2 in the present case) and 
which is located at a time-dependent radius 7v(t). If 
the radial velocity is small with respect to the azimuthal 
velocity (as is true during most of the plunge) we can 
think of a source with an adiabatically varying frequency 
mTt{t) located at the adiabatically varying radius r*(t). 
Replacing this behavior in equation (|10p allows one to 
approximately replace the repeated time derivative by 
the square of the instantaneous source frequency, i.e. 
by a factor (— imfl(t)) 2 = —(mfl(t)) 2 . This yields a 
Schroedinger-type (radial) equation with effective poten- 
tial V^'°\n) ~ {mVL{t)f. In the upper panel of Fig. gj 
we compare the evolution, as a function of dynamical 
time t, of the square roots of the two contributions to 

this effective radial potential: \Jv^°\r^{t)) and mfi(i) 
( for the case t = m = 2). We see that the (mil(t)) 2 
term is always significantly smaller than the V^°' (r* (t)) 
one. This means that during all the plunge (including at 
the end, when £l(t) reaches its maximum value), we can 
think (in first order) in terms of a near-zone approxi- 
mation, where we can neglect (2f2(i)) 2 and compute the 
asymptotic waveform by convolving the source by a so- 
lution of the static, radial Reggc- Wheeler- Zcrilli equa- 
tion. One can check that this leads essentially to the 
usual (Schwarzschild-coordinate) quadrupole formula for 
the waveform. In this approximation there is strictly no 
excitation of the QNM modes by the adiabatically vary- 



ing source. One therefore concludes (following the line 
of reasoning of Ref. [24[ where it was explained how a 
certain type of oscillating integral which exponentially 
vanishes in the adiabatic approximation acquired a non- 
zero value due to a localized non-adiabatic behavior of 
the integrand) that the QNM modes will be excited es- 
sentially only around the maximum of the exciting fre- 
quency mQ(t), because this is the moment when the 
source excitation is non-adiabatic. Looking now at the 
bottom panel of Fig. g[ we see the trajectory of the ex- 
citing frequency on the positive real axis, from zero up 
to the maximum 2J7 max = 0.2703. The black dots in 
the complex frequency plane represent the first QNM 
modes of a Schwarzschild black hole. One understands 
intuitively (and this can be put in equations similar to 
those derived in [24|]) that the amount of excitation, by 
the non-adiabatic behavior of the excitation frequency 
near its maximum, of the various QNM modes will pri- 
marily depend on their "distance" to the critical (real) 
exciting frequency +2Sl max = +0.2703. This reasoning 
indicates that the QNM modes which will be primarily 
excited are the first few positive frequency modes, i.e. 
those appearing on the right part of the bottom panel of 
Fig. g[ Summarizing in slightly different physical terms, 
we can think of the black hole as a resonating object hav- 
ing the resonance spectrum showed in Fig. |4j When the 
source varies adiabatically, so that it corresponds to a 
nearly pure frequency +mtt which is much smaller than 
any of the resonance frequencies it excites only negli- 
gibly the resonance modes. However, when the source 
varies non-adiabatically (near its maximum frequency), 
its Fourier transform contains a halo of frequencies (due 
to non-adiabaticity) centered around the instantaneous 
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FIG. 4: Illustration of some of the important physical features 
of the excitation of quasi-normal modes by a small-mass-ratio 
coalescing binary system. The upper panel shows that, dur- 
ing all the plunge, one remains in a "quasi-stationary" near- 
zone regime (which does not excite QNMs). The crucial fea- 
ture which can excite the QNMs is the non-adiabaticity of the 
evolution of the exciting orbital frequency near its maximum 
(similarly to Ref. [241]'). The lower panel illustrates why the 
(rather short) non-adiabatic behavior of 2Q(t) near its maxi- 
mum preferentially excites the positive-frequency QNMs. 



frequency, and this halo can now "extend the reach" of 
+mf2 max and thereby excite the black hole resonances 
which lie within this "halo" (i.e. the first few QNM 
modes whose frequency has the same sign as +mfi max ). 



IV. EOB-TYPE APPROXIMATE WAVEFORMS 

A. Methodology: a matched waveform 

Let us now turn to the construction of complete EOB- 
type gravitational waveforms obtained by matching an 
analytical inspiral-plus-plunge waveform to a ring-down 
one. For the reasons explained above, the crossing of the 
light ring (actually the non adiabatic behavior around 
the moment when Q(t) reaches it maximum) corresponds 
to an abrupt triggering of the QNMs of the black hole 
(this was realized long ago [37] for the case of the ra- 
dial plunge) . The EOB approach takes into account this 
abrupt change in the underlying physics by an abrupt 
change in the analytical description of the system: be- 
fore the crossing of the light ring one describes the bi- 
nary system as two point masses with EOB relative dy- 
namics, while after the crossing one replaces the binary 
system by a single distorted black hole (as in the "close- 
limit" approximation of colliding black holes Q). As 
we were mentioning in the introduction, at the practical 



level this dual description implies that one should match 
together (around the light ring) a plunge waveform ob- 
tained with a quadrupole-type formula to a superposition 
of black hole QNMs. This general philosophy has already 
been implemented, with increasing sophistication in Refs. 
[E [Hi [l7], [U, H3| • Our aim here is to further improve the 
implementation of the EOB approach by studying in de- 
tail each one of its building blocks, and by comparing 
the resulting analytical waveform to the exact numerical 
waveforms discussed above. 



B. An improved, resummed 3PN-accurate 
inspiral-plus-plunge waveform. 

We start by introducing here a new, improved EOB- 
type inspiral-plus-plunge waveform that goes beyond pre- 
vious work in two directions: (i) it includes higher post- 
Newtonian (PN) effects which are resummed in a novel 
way (factorized formulation, exact resummation of tail 
effects, 3PN- resummed non-linear PN effects), and (ii) 
it includes non-quasi-circular (NQC) corrections that are 
expected to be relevant below the LSO. Details of our 
construction will be explained separately 36]. Up until 
the light ring (r > 3), we use a (£,m) — (2,2) gravita- 
tional waveform defined by means of the following EOB- 
type improved quadrupole formula 



Pr, 



22 



. u Pr, 
10 

rQ 



-2i<p 



(14) 



Here, F 22 includes four types of PN effects, while 
the braces represent NQC corrections. The standard 
quadrupole formula (for a point-particle source in quasi- 
circular motion) corresponds to the case where F 22 = 1 
and a = b = (see e.g. [IH; the numerical coefficient 
\J 7r/30 corresponds to our using the Zerilli-Moncrief even 

quantity \& 2 2 )■ 

The NQC factor added above (within braces) is more 
general than the general time-even NQC factor consid- 
ered in (fTTj) because we describe here non time-even NQC 
effects (as can easily be seen on the usual quadrupole for- 
mula written with two time-derivatives). The numerical 
coefficients a and b thus parametrize non-quasi-circular 
corrections to the modulus and to the phase of the wave. 
They will be fixed below by comparison to the "exact" 
Zerilli-Moncrief waveform. The novel F22 PN-improving 
factor is given as the product of four terms 

F 22 (t) = #(f)T(2, 2Sl(t))e i5 ^f 22 (x(t)) . (15) 

This factorization of PN effects is, in itself, a new way 
of partially resumming PN effects (by contrast to usual 
representations of PN effects which mix the expansion 
of the four separate factors above). The first factor is 
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the "conserved" relativistic energy of the (effective) mov- 
ing source. The second factor T(2, 20) is the partic- 
ularization to I = m = 2 of a function T(£, k) (with 
Mk = k — mO) that represents a "tail" correction, 
namely the contribution due to backscattering of (multi- 
polar) waves off the background curvature. This factor 
is the exact resummation of an infinite number of "lead- 
ing logarithms" , and is derived, similarly to the work in 
Ref. [39j, by using Coulomb wave functions. It is explic- 
itly given by 



C. Ring-down waveform 

After the crossing of the light ring, the EOB approach 
consists in replacing the two-body system by a single dis- 
torted black hole. At the level of waveforms, this means 
that one should match, around the light ring, the plunge 
waveform Eq. (|14[) to a ring- down one made of a super- 
position of (I = 2, m = 2) black hole QNMs, say 



^ringdown^ 



E 



C+e 



c- 



(18) 



T(£,k) 



r(* + l-2ifc) „ 



T(t + 1) 



e -k e 2ik log(2fc& scal<i ) 



(16) 



where we shall take & S caic = 2M. e 1 * 22 is an additional 
phase correction, with S22 = 7fi/3 to lowest PN or- 
der. The final factor f22(x) represents the remaining 
(essentially non-linear) PN effects. Using Ref. [4(| it 
can be computed in the test-mass limit, during the inspi- 
ral, as a Taylor expansion in the PN-ordering parameter 
x = l/r = ft 2 / 3 . At 3PN 8 , it reads 



f T^or {x) = 1 __ i _ 



43 
21 : 
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1712 
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log(2) 



21428357 
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428 



105 



log(x) 



(17) 



where 7 = 0.577216 ... is the Euler constant. To improve 
its convergence properties in the strong-field-fast-motion 
regime, expression (|17p is further Pade resummed in a 
proper way. 

When deriving the factor /22O1O during the inspiral 
phase, the argument x that it depends upon is "degen- 
erate" in the sense that it can be expressed in various 
ways in terms of the dynamical variables of the system, 
namely: x = u 2 = {rfl) 2 = ft 2 / 3 = l/r. However, these 
various expressions start differing from each other dur- 
ing the plunge. This opens the issue of choosing the 
"best" dynamical expression for the PN-ordering argu- 
ment x during the plunge, in the sense of giving the best 
approximation to the "exact" waveform (minimizing the 
need for NQC corrections). We have found that choosing 
x = l/r (the Newtonian gravitational potential) yields 
good results, and this is the choice we shall make here. 
We postpone to a further publication a thorough analy- 
sis of the effect of different x arguments, as well as the 
contrast between the Pade-resummed (that we shall use 
here) and Taylor-expanded versions of f%2{ x )- 



where — ot n ± iw„ are the positive/negative complex 
QNM frequencies and are the corresponding ampli- 
tudes. Here, a n and co n indicate the damping time and 
the oscillation frequency of each mode respectively, and 
n = 0, 1, 2, N — 1 labels the overtone number (n = 
denoting the fundamental mode). 

In this paper we propose an improved treatment of the 
definition of the ring-down waveform Eq. (|18[) . as well as 
of the matching to the plunge waveform Eq. (fT4")l . First, 
in view of both the numerical results and the theoretical 
arguments presented above, we know that the contribu- 
tion of the positive-frequency modes to the waveform is 
much larger than that of the negative-frequency ones. 
Therefore, we can set C~ = in first approximation. 
As a result, the instantaneous gravitational wave fre- 
quency Mwf 2 OB = -3(*f 2 OB /*f 2 OB ) of our "matched" 
ring-down waveform will grow until it reaches a "flat" 
saturation region, without oscillations. Though this is 
qualitatively slightly different from the behavior of the 
exact instantaneous frequency (see Fig. [3]) we will see 
that this gives an excellent approximation to the phasing 
(and actually a better one that if one matches to both 
positive and negative modes in a democratic way). 

Our problem is thus to determine the coefficients C+ 
from the knowledge of ^2 2 ungc for r > 3. Here again 
we introduce some significant improvements on previ- 
ous implementations of the EOB philosophy, based on 
the knowledge acquired from detailed comparisons with 
our "exact" small-mass-ratio numerical results. A triv- 
ial improvement will consist in using a larger number N 
of QNMs (including the fundamental mode) than pre- 
vious work. Instead of N = 1 [H, |3j, N = 2 pi], or 
N = 3 [13, HU, we shall use here N = 5 (positive- 
frequency) quasi-normal modes (including the fundamen- 
tal mode) [3J,l35j . Our studies showed that N = 5 gives a 
significantly better fit than N = 4 or less, and that larger 
values of N yield only a rather marginal improvement. 



D. Determining the next-to-quasi-circular (NQC) 
corrections to the plunge waveform. 



We note that for i/ < 1 this series is known to higher PN orders. 
We truncate it here to the 3PN level which is the highest order 
at which we could (and did) derive it in the comparable-mass 
case. We leave to a future publication a detailed discussion of 
the derivation of these results. 



First of all 'J'^" 118 ' 6 ! Eq. (|14j) . has two free parameters a 
and b that should be specified so to minimize the differ- 
ence with the real signal. The presence of two separate 
NQC real parameters is useful because the waveform is a 
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complex number, and we wish to best fit both its modu- 
lus and its phase 9 . The modulus of the numerical wave- 
form reaches a maximum at the (shifted) retarded time 
u/(2M) — 299.6 (i.e. just before the maximum of the 
orbital frequency, or the crossing of the light ring, which 
both occur near u/(2M) — 300.9). We can then con- 
strain the NQC parameters a and b so as to ensure that 
the maximum of the modulus of the analytical waveform 
([14]) sits on top of that of the exact modulus. Actually, 
we can even do this analytically. Using the quasi-geodesic 
approximate description of the plunge, and imposing that 
the modulus of the analytical waveform (|14[) is maximum 
at the light-ring gives, in first NQC approximation, the 
constraint 2a + b 2 = 1. We then get a second constraint 
(mainly on the coefficient b) by requiring that the slope 
of the instantaneous analytical frequency Mcu^'^it) be 
as close as possible to the exact one Muj22(t) around the 
light ring. The two above requirements fix the values of 
the NQC parameters to be: a = 0.438 and b = 0.355. It 
is remarkable that these two numerical values (numeri- 
cally determined, by trial and error, as best eye-fitting 
values) do satisfy quite accurately the (first-order) ana- 
lytically determined constraint mentioned above: indeed, 
2a num + 6„ um = 1.002. This shows the power of the EOB 
approach for building a purely analytical self- consistent 
matched waveform, even in absence of knowledge about 
the exact signal. 



E. Matching the inspiral-plus-plunge waveform to 
the ring-down one on a "comb". 

The last and most crucial new ingredient introduced 
here consists in matching the plunge and ring-down wave- 
forms on a finite time interval [t rn — 8/2, t m + 8/2], in- 
stead of simply imposing the continuity of the two func- 
tions, and a certain number of their derivatives, at one 
moment t matc h (as done in [J, Qj], [U, H3]). We found 
that this matching on a finite interval (chosen in the 
immediate vicinity of light-ring crossing) was quite ef- 
fective for accurately reproducing the smooth rise of 
the instantaneous GW frequency during the whole tran- 
sition. More precisely, a good choice for the match- 
ing interval [t m — 8/2, t m + 8/2] consists in putting its 

center t m at the maximum of the modulus of \&22 UnSC j 
i.e. at t/(2M) = 299.6, and choosing his full length 
8 — At as S/M ~ 7.2. The matching on the interval 
[t m — 8/2, t m + 8/2] is most easily done by discretiz- 
ing this matching interval into N (here N = 5) equally 



9 In view of this, a useful alternative to the a, b NQC 
factor in braces used above, is to use instead a factor- 
ized complex NQC waveform correction of the form, (1 + 
a 'Pr /( r ^) 2 )) exp(+i&'p r , /(rf2)), in which a' affects only the 
modulus, and b' only the phase. The (approximate) analytic 
constraint on the modulus discussed below then becomes simply 
2a' = 1. 



spaced points, [ti, t<z, . . . , ijv] (with t\ — t m —8/2 and 
ijy = t m + 8/2) and requiring that the two complex func- 
tions (fl4|) and (fT8|) coincide at these N points. This 
gives N (linear) complex equations for the N complex 
unknowns C+. In other words we are using a 'matching 
comb'. 



F. Comparing analytical and numerical waveforms 

Finally, we define our complete EOB matched wave- 
form (from inspiral to ring-down) as 



EOB/,a _ 
22 



(t) = 0(t r 



(19) 



where 8(t) denotes Heaviside's step function 10 . 

The left panel of Fig. [3] compares the (real and imagi- 
nary parts of the) matched EOB-type waveform (fl"9|) (red 
line) obtained by means of the matching procedure ex- 
plained above, to the "exact" numerical one (black line). 
The agreement between the two is clearly excellent. Note, 
in particular, how the inclusion of the novel 3PN-accurate 
correction factor F22 allows the analytical inspiral wave- 
form to coincide in phase and in amplitude with the exact 
signal. As we shall illustrate below, Newtonian-accurate 
"restricted" waveforms do a much poorer job: though 
they describe the phasing rather accurately, they over- 
estimate (already during the inspiral) the amplitude by 
a significant factor. To investigate further the quality of 
this agreement, we also plot the modulus of the wave- 
forms and the instantaneous frequency (right panel of 
Fig. [3]) . One sees again on this figure the very accurate 
way in which the improved EOB signal reproduces both 
the frequency and the modulus of the exact signal dur- 
ing the inspiral and most of the plunge (which starts at 
u/2M ~ 240). It is only at the end of the plunge, just be- 
fore light-ring crossing, (and during ring-down) that one 
notices a difference (smaller than about 5%) in modulus. 
To further study the difference in frequency of the two 
signals, we compare in Fig. [5] the time evolution of the 
phases of the two signals (with the convention W = e~ 1 ^ , 
so that (f> = +lu): the matched instantaneous gravita- 
tional wave phase </>f 2 OB versus the exact one 4^22^' '■ 
These quantities are computed from time integration 
of the instantaneous frequencies from an initial time to 
starting with </>o = 2ipo where ipo is the initial orbital 
phase. The visual agreement between the exact phase 
(black line) and the matched phase (red line) is so good 
that the two lines can barely be distinguished. At a quanti- 
tative level, the phase difference A^f^ 



EOB _ i,( cxact ) 
22 ^22 



10 If one wanted to have a C°° transition between the two wave- 
forms one could replace 8(t — t m ) by one of Laurent Schwartz's 
well-known smoothed step functions (or "partitions of unity") 

e c ((t-t m )/s) 
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FIG. 5: Analytic and numerical gravitational wave phases 
versus (shifted) retarded time. As they are very nearly super- 
posed, the inset shows the difference between the two phases 
(computed, as a check, by means of two independent meth- 
ods). 



FIG. 6: Comparison between the moduli and the phases 
of two waveforms: the "exact" numerical one, and a non- 
optimized analytic one based on a coarser implementation of 
the EOB philosophy (no PN corrections, no explicit NQC cor- 
rections, N = 3 QNMs, nearly instantaneous matching at the 
light ring). See text for details. 



is an important diagnostics which is shown in the in- 
set of of Fig. [5] (after adding a constant offset). This 
has been computed by two separate methods. The first 
one (red dashed line) computes A^f™ by straightforward 
time integration of the difference between Mwf 2 OB an d 
Mu>22- The second computation (black solid line) does 
not involve any time integration, it just evaluates the ar- 
gument of the ratio of the two complex numbers v Pf 2 ° B 



and * 



(exact) 
22 



^f 2 OB |^ 2 ° 2 XaCt) | 
^(cxact) |vjjEOB| 



(20) 



As shown in the inset of Fig. [5] there is an excellent con- 
sistency between the two methods. 

Fig. [5] exhibits our most important final result: the 
maximum phase difference, during the whole transition 
from inspiral to plunge, between our improved EOB sig- 
nal and the exact one is of the order of ±1.1% of a cycle. 
Note that the ripples in the late-time part of At/if^ are 
(obviously) related to the absence of negative frequency 
modes in the matched waveform. 

This excellent analytical description of the exact wave- 
form (phasing to ±1.1% of a cycle, and a rather accurate 
modulus, see Fig. [3]) was obtained by combining the sev- 
eral improvements to the EOB method introduced above, 
and by using the specific numerical values of the new pa- 
rameters aim, a, 6, t m , S quoted above. Note, however, 
that there remains some flexibility in the numerical val- 
ues of these parameters. We leave to future work a de- 
tailed discussion of this issue. Let us only note here that 
the "good choices" of these parameters (minimizing the 



maximum (absolute) dephasing n ) are correlated, and 
that the values we have used above are near optimal in 
the following sense: if, for example, we keep fixed the 
matching comb parameters t m and 5, and vary a in the 
interval [0.3, 0.8], the best possible complementary values 
of b that we could find (varying in the interval [0.2, 0.5]) 
were found to lead to minimax dephasings larger than 
1.1% , but remaining smaller than 3%. We also explored 
the sensitivity to the choice of the matching comb. For 
instance, if we keep fixed only the center of the matching 
interval at t m / (2M) — 299.6, and increase S/M by a fac- 
tor 2 (to S/M = 14.4) and then look for the best possible 
values of the parameters a and b, we find that the min- 
imax dephasings are increased from 1.1% to about 2%. 
On the other hand, if (still keeping fixed t m ) we decrease 
S/M to 1.8, and then look for the best possible a and b, 
we find that the minimax dephasings are increased from 
1.1% to 2.5%. Similar sensitivities are also reached if we 
keep fixed S/M (to the near optimal value 7.2), but move 
the center of the interval: for example, t m /(2M) = 298 
leads to a minimax of 2%, while t m /(2M) — 301 gives 
4%, so that we obtain correspondingly less accurate de- 
criptions of the waveform. 

On the other hand, as we shall discuss in detail in a 
separate publication, even much simplified versions of our 
matching technique can lead to rather impressively good 
results, notably for the phasing. To quote one particu- 



11 We shall refer to the minimum value of the maximum absolute 
value of the dephasing (over the whole inspiral-plunge-merger 
period) as the "minimax" dephasing. 
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larly simple example, the standard quadrupole formula in 
the quasi-circular approximation (in the sense of Eq. (|14[) 
with F22 = 1, a = b = and N = 5), still gives a quite 
tiny dephasing, A0|^/(27r) < 2%, though the evolution 
of the modulus is less well reproduced. 

Finally, we pedagogically illustrate in Fig. [6j in a "con- 
trasting" manner, the effect of the various improvements 
introduced above by showing how the waveform compar- 
ison of Fig. [3] gets modified when not using them, but 
using instead the type of coarser implementation of the 
EOB philosophy used in previous work [J, i, EE El El . 
More precisely, the EOB-type analytical waveform used 
in Fig. [6] was obtained by: (i) using (as originally pro- 
posed in Ref. 1]) the following (Newtonian-order and 
Kepler-law-assuming) restricted quadrupole waveform 



lUPUp, El- The current 
of this method have 



^ 2 K (i) = -4 M 



30 



(21) 



without any explicit PN (F22) corrections, nor any NQC 
(a, b) corrections; (ii) only 3 (positive-frequency) QNMs; 
and, rather crucially, (iii) by matching the plunge and 
ring-down waveforms in a very small interval (5/M = 0.9 
instead of our preferred 7.2) around the maximum of the 
orbital frequency. [Indeed, the matching of the two wave- 
forms and their derivatives at a sharply defined moment 
is equivalent to considering the <5 > limit of our comb- 
matching technique] . By contrasting Fig. [5] with the right 
panel of Fig. [3] we see that: (i) the modulus of the ana- 
lytical waveform is now distinctly larger than the exact 
one during the inspiral (because of the lack of PN correc- 
tions) ; (ii) the modulus becomes significantly larger than 
the exact one at the end of the plunge (because of the 
use of the Kepler-law-assuming cx Q 2 / 3 , which, as pointed 
out in |24[ tends to overestimate the amplitude) ; (iii) the 
post-matching analytical frequency jumps up from the 
(doubled) maximum orbital frequency significantly more 
vertically than before, thereby decoupling too soon from 
the exact frequency, and accruing a larger dephasing than 
before (because of the too localized matching, and - to a 
lesser degree - the use of only 3 QNMs). Note, however 
that (as exhibited in the inset) this "coarser" EOB-type 
implementation still succeeds in following the phase of 
the exact signal to a rather impressive ±3.5% of a cycle. 



V. CONCLUSIONS 

We have presented a first attack on the problem 
of computing analytic template waveforms, describing 
the complete binary black hole coalescence phenomenon 
(from early inspiral to late ring-down) that are both ef- 
fectual and faithful (in the terminology of Ref. [20] ) . We 
have restricted our analysis to the small-mass-ratio limit 
without spin, a problem that is amenable of a semi- 
numerical investigation based on a synergic use of black 




hole perturbation theory and PN-based EOB results [2£ 

Our analytic approach is based on the Effective One 
Body (EOB) method 
implementations [l], 0. 
been recently shown [17j to define effectual templates 
for the coalescence of comparable-mass (non-spinning) 
systems in the full mass range of interest for ground- 
based GW detectors. Our aim is to improve the imple- 
mentation of the EOB method so as to construct high- 
accuracy, faithful templates. Our strategy towards this 
aim is to understand, element by element, the physics be- 
hind each feature of the waveform, by comparing general- 
ized, multi-parameter EOB-based waveforms to "exact" 
waveforms obtained by numerical relativity approaches. 
In this first paper, we have tapped the information con- 
tained in small-mass-ratio waveforms, which can rather 
easily be computed very fast and with high-accuracy. 

Our main achievements are: (i) a better description 
of radiation reaction during the plunge (by taking into 
account non-quasi-circular effects); (ii) a refined analyt- 
ical expression for the (£, m) = (2,2) waveform (which 
includes a resummation of tail effects, a 3PN-accurate 
"non-linear" amplitude correction, and the inclusion of 
non-quasi-circular corrections); (iii) an improved tech- 
nique for matching the plunge and ring-down waveforms 
(which is based on an improved understanding of the 
physical elements entering this matching: origin of the 
excitation of quasi- normal modes, sign of their frequen- 
cies, use of an extended matching interval). The imple- 
mentation of these improvements has led us to the ex- 
plicit computation of EOB-type analytic waveforms that 
accurately reproduce (both in its phase and its ampli- 
tude) the exact waveform during the whole transition. 
In particular, the phase difference stays within ±1.1% of 
a cycle. Such a good result in the small mass-ratio case 
rises the hope that a comparably good agreement could 
also be achieved in the comparable-mass case. This will 
be the object of future work [4lj . 
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